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Abstract 

A partial wave analysis of the centrally produced 7r + 7r~ channel has been performed in pp 
collisions using an incident beam momentum of 450 GeV/c. An unambiguous physical solution 
has been found. Evidence is found for the / (980), / (1300), / (1500) and /j(1710) with J = 
in the the S-wave. The p(770) is observed dominantly in the P^-wave and the /2(1270) is 
observed dominantly in the .D^-wave. In addition, there is evidence for a broad enhancement 
in the D-wave below 1 GeV. 
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Lattice gauge calculations indicate that the lowest lying glueball should be a scalar and 
be in the mass range 1500-1700 MeV [l|]. However, after more than 20 years of searching 
for glueballs there is emerging evidence that a pure scalar glueball may never be observed 
due to the fact that its finite width means that it mixes with nearby isoscalar qq states with 
J PC = ++ 0. This scenario gives rise to several states whose properties are complicated and 
lie somewhere between conventional qq states and glueballs. Hence a systematic study of all 
isoscalar J PC = ++ states is required. The nn system has long been a profitable place to study 
jGjpc _ q+0 ++ states and to date very high statistics studies have been performed using 7r 
induced interactions || and pp annihilations ||. This paper presents a study of the -k + tt~ final 
state formed in central pp interactions which are predicted to be a source of gluonic final states 
via double Pomeron exchange ||. 

The reaction 

pp -»• Pf(ir + ir~)p s (1) 

has been studied at 450 GeV/c. The subscripts / and s indicate the fastest and slowest 
particles in the laboratory respectively. The WA102 experiment has been performed using the 
CERN Omega Spectrometer, the layout of which is described in ref. ||. Reaction (1) has 
been isolated from the sample of events having four outgoing charged tracks, by first imposing 
the following cuts on the components of the missing momentum: |missing P x \ < 14.0 GeV/c, 
|missing P y \ < 0.16 GeV/c and |missing P z \ < 0.08 GeV/c, where the x axis is along the 
beam direction. A correlation between pulse-height and momentum obtained from a system of 
scintillation counters was used to ensure that the slow particle was a proton. 

The method of Ehrlich et al. J7|, has been used to compute the mass squared of the two 
centrally produced particles assuming them to have equal mass. The resulting distribution is 
shown in fig. |l|a) where a clear peak can be seen at the pion mass squared. A cut on the Ehrlich 
mass squared of —0.3 < M\ < 0.2 GeV 2 has been used to select a sample of 5.15 million 
events. 

Fig. |I]b) shows the Pf^ + effective mass spectrum where a clear A ++ (1232) can be observed. 
A smaller A°(1232) signal can be seen in the p/vr~ effective mass spectrum shown in fig. 0c). 
The A(1232) signal has been removed by requiring M(j)fir) > 1.5 GeV. Due to the trigger 
requirements there is little evidence for A production in the p s 7i effective mass spectra (not 
shown). However, for symmetry purposes a cut of M(p s 7i) > 1.5 GeV has also been applied. 

As can be seen from fig. |l]b) there is also some evidence for higher mass A ++ or iV* produc- 
tion. In order to investigate the effect of this on the following results an analysis has also been 
performed requiring that M(piv) > 2.0 GeV. It is found that the higher mass proton excitations 
do not significantly influence the results. 

The Feynman xp distributions for the slow particle, the 7r + 7r" system and the fast particle 
are shown in fig. |l]d). As can be seen the 7r + 7r~ system lies within \xp\ < 0.25. 

The resulting centrally produced n + n~ effective mass distribution is shown in fig. |l|e) and 
consists of 2.87 million events. As can be seen there is evidence for a small p°(770) signal, some 
/2(1270) and a sharp drop at 1 GeV which has been interpreted as being due to the interference 
of the /o(980) with the S-wave background ||. 

A Partial Wave Analysis (PWA) of the centrally produced 7r + 7r~ system has been performed 
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assuming the 7r + 7r~ system is produced by the collision of two particles (referred to as exchanged 
particles) emitted by the scattered protons. The z axis is defined by the momentum vector of 
the exchanged particle with the greatest four-momentum transferred in the 7r + 7r~ centre of 
mass. The y axis is defined by the cross product of the two exchanged particles in the pp centre 
of mass. The two variables needed to specify the decay process were taken as the polar and 
azimuthal angles (9, <p) of the n~ in the centre of mass relative to the coordinate system 

described above. 

The acceptance corrected moments, defined by 

J(ft) = £t io y L °(ft) + 2 J2 t LM Re{Y^(n)} (2) 

L L,M>0 

have been rescaled to the total number of observed events and are shown in fig. |2|. As can be 
seen the moments with M > 2 and L > 4 are small (i.e. i 43 , £ 44 , £ 50 and £ 60 ) and hence only S, 
P, and D waves with m < 1 have been included in the PWA. 

An interesting feature of the moments is the presence of a structure in the L = 4 moments 
for ir + 7T~ masses below 1 GeV which indicates the presence of D-wave. This type of structure 
has not been observed in n induced reactions ||. In order to see if this effect is due to ac- 
ceptance problems or problems due to non-central events, we have reanalysed the data using a 
series of different cuts. Firstly, we required that M(p7r) >2.0 GeV; however, after acceptance 
correction the moments were compatible with the set for M(jm) >1.5 GeV showing that diffrac- 
tive resonances in the range 1.5 < M(pn) < 2.0 GeV have a negligible effect on the moments. 
We have also required that the rapidity gap between any proton and n in the event is greater 
than 2 units. Again the resulting acceptance corrected moments do not change. In order to 
investigate any systematic effects we have also analysed the central 7r°7r° data and a similar 
structure is also found in the L = 4 moments §|. 

It is interesting to compare the above result with previously published data on other centrally 
produced ttty systems. In preliminary results from the E690 experiment at Fermilab activity 
at a similar level is observed in the £ 4 o moment below 1 GeV llOfl . There is also evidence for 
this structure in the 7r°7r° data from the NA12/2 experiment [IT]. The AFS experiment at the 



CERN ISR also observed that the £ 4 o moments deviated from zero in this mass region; however, 
in their analysis they claimed this deviation was due to problems of the Monte Carlo simulating 
low energy tracks ]12| . 



This structure does indeed seem to be a real effect which is present in centrally produced 
7T7T systems. It has recently been suggested |L3], [L4]] that central production may be due to the 
fusion of two vector particles and this may explain why higher angular momentum systems can 
be produced at lower masses. 



The amplitudes used for the PWA are defined in the reflectivity basis |I5] . In this basis the 
angular distribution is given by a sum of two non-interfering terms corresponding to negative 
and positive values of reflectivity. The waves used were of the form with J = S, P and D, 
m = 0, 1 and reflectivity e = ±1. The expressions relating the moments (£lm) and the waves 
(Jm) are given in table |TJ. Since the overall phase for each reflectivity is indeterminate, one wave 
in each reflectivity can be set to be real (Sq and Pj 4 " for example) and hence two phases can be 
set to zero {<f> s - and (p p + have been chosen). This results in 12 parameters to be determined 
from the fit to the angular distributions. 
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The PWA has been performed independently in 20 MeV intervals of the 7r + n mass spec- 
trum. In each mass an event-by-event maximum likelihood method has been used. The function 

N 

F = -Y, ln{I(n)} + £ t LM e LM (3) 

i=l L,M 

has been minimised, where N is the number of events in a given mass bin, €lm are the effi- 
ciency corrections calculated in the centre of the bin and tiM are the moments of the angular 
distribution. The moments calculated from the partial amplitudes are shown superimposed on 
the experimental moments in fig |2|. As can be seen the results of the fit well reproduce the 
experimental moments. 

The system of equations which express the moments via the partial wave amplitudes is 
non-linear that leads to inherent ambiguities. For a system with S, P and D waves there are 
eight solutions for each mass bin. In each mass bin one of these solutions is found from the 
fit to the experimental angular distributions while the other seven can then be calculated by 



the method described in ref. |15[]. In order to link the solutions in adjacent mass bins, the 
real and imaginary parts of the Barrelet function roots are required to be step-wise continuous 
and have finite derivatives as a function of mass [[[(J. By definition, all the solutions give 
identical moments and identical values of the likelihood. The only way to differentiate between 
the solutions, if different, is to apply some external physical test, such as requiring that at 
threshold that the S-wave is the dominant wave. 

The four complex roots, Zj, after the linking procedure are shown in fig. |3|. As can be seen 
the real parts are well separated and hence it is possible to identify unambiguously all the 
PWA solutions in the whole mass range. In addition, the zeros do not cross the real axis and 
hence there is no problem with bifurcation of the solutions. Near threshold the P-wave is the 
dominant contribution for five solutions, another one is dominated by D-wave and another has 
the same amount of S-wave and P-wave. These seven solutions have been ruled out because 
the 7r + 7r~ cross section near threshold has been assumed to be dominated by S-wave. The 
remaining solution is shown in fig. |j. 

The S-wave spectrum shows a clear threshold enhancement followed by a sharp drop at 
1 GeV. There is clear evidence for the p(770) in the Pq wave and for the /2(1270) in the Dq 
wave. 

In order to obtain a satisfactory fit to the Sq wave from threshold to 2 GeV it has been 
found to be necessary to use three interfering Breit-Wigners to describe the / (980), / (1300) 
and /o(1500) and a background of the form a(m — m t h) b exp(—cm — dm 2 ), where m is the 
7r + 7r~ mass, mth is the n + n^ threshold mass and a, b, c, d are fit parameters. The Breit- 
Wigners have been convoluted with a Gaussian to account for the experimental mass resolution 
(o" = 4 MeV at threshold rising to 22 MeV at 2 GeV). The fit is shown in fig. |3]c) for the entire 
mass range and in fig. |3|d) for masses above 1 GeV. The resulting parameters are 

/o(980) M = 982 ± 3 MeV, r = 80 ±10 MeV 
/o(1300) M = 1308 ± 10 MeV, T = 222 ± 20 MeV 
/o(1500) M = 1502 ± 10 MeV, T = 131 ± 15 MeV 

which are consistent with the PDG |T7J values for these resonances. As can be seen, the fit 
describes the data well for masses below 1 GeV. It was not possible to describe the data above 
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1 GeV without the addition of both the / (1300) and / (1500) resonances. However, even with 
this fit using three Breit-Wigners it can be seen that the fit does not describe well the 1.7 GeV 
region. This could be due to a n + Ti~ decay mode of the fj{YJVS) with J = 0. Including a 
fourth Breit-Wigner in this mass region decreases the y 2 from 256 to 203 and yields 

/j(1710) M = 1750 ± 20 MeV, T = 160 ± 30 MeV. 



parameters which are consistent with the PDG [17] values for the /j(1710). The fit is shown 
in fig. [5]e) for masses above 1 GeV. 

In conclusion, a partial wave analysis of a high statistics sample of centrally produced 
7r + 7r~ events has been performed. An unambiguous physical solution has been found. The 
S-wave is found to dominate the mass spectrum and is composed of a broad enhancement 
at threshold, a sharp drop at 1 GeV due to the interference between the /q(980) and the S- 
wave background, the / (1300), the / (1500) and the /j(1710) with J = 0. The p(770) is 
observed dominantly in the P ~-wave. The D-wave shows evidence for the /2(1270) and a broad 
enhancement below 1 GeV. It is interesting to note that the /2(1270) is produced dominantly 
with m = 0. There is no evidence for any significant structure in the D-wave above the /2(1270) 
i.e. no evidence for a J = 2 component of the /j(1710). 
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Table 1: The moments of the angular distribution expressed in terms of the partial waves. 



v/i^oo = \So\ 2 + \Po\ 2 + \p^\ 2 + \p^\ 2 + \Do\ 2 + \D^\ 2 + \Dt\ 2 

VAnho = 2|S -||Po1 cos Os- - 4> P - ) + 7g I ?o 1 1 D o I cos {4> P - ~ 4> D ~ ) 

+ 2 -${\Pr\\Di\cos(<f> Pi - - 4> D -) + \P+\\Dt\cos{<t> Pt - </> D +)} 

v^Ftii = v^|5 , o -||Pr|c OS (0 S(7 -0 Pr )-^|Pfi|D o -|c OS (0 Pr - ( /) D -) 
+^l^rlPr|aw(0p--^ r ) 

+ 2 1^0 II A) I COS(0 5 - - 4> D -) 

+v / 2|5' o -|| J Dr|cos(0 5 - -0 D -) 

= ^(IPfP-IA+n + ffdPrl 2 -!^! 2 ) 

= -^{|pri|p>r|co S (0 Pi --0 Dr ) + |p 1 + ||D 1 -|co S (0 Pi+ -0 D +)} 

+ ^I^o1IA1co S (0 p --0 d -) 

= Jg|Prl|Po1co S (0 Pr -0 D -) + ^|p o -|pr|co S (0 P --0 Di -) 

^32 = ^{iPfllPrlcO^^- - Dr ) - |P+||P>+| COS (0 p + - 

= f|p l 2 -f(l^rl 2 + |Pi-| 2 ) 
v/i^ = ^f(|prl 2 -l^l 2 ) 



7 



Figures 

Figure 1: a) The Ehrlich mass squared distribution, b) the M (pfTc + ) and c) the M(p/7r _ ) mass 
spectra, d) The xp distribution for the slow particle, the 7r + 7r~ system and the fast particle, 
e) The centrally produced ti + ti~ effective mass spectrum. 



Figure 2: The ^/4nt LM moments from the data. Superimposed as a solid histogram are the 
resulting moments calculated from the PWA of the ti + ti~ final state. 



Figure 3: a) The Real and b) Imaginary parts of the roots (see text) as a function of mass 
obtained from the PWA. c) and d) The 7r + 7r - Sq wave with fit described in the text using three 
Breit-Wigners. e) The 7t + tt~ Sq wave with fit described in the text using four Breit-Wigners. 



Figure 4: The physical solution from the PWA of the tt + tt final state. 
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